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We use angular size versus redshift data for galaxy clusters from Bonamente et al. [l[ to place 
constraints on model parameters of constant and time-evolving dark energy cosmological models. 
These constraints are compatible with those from other recent data, but are not very restrictive. A 
joint analysis of the angular size data with more restrictive baryon acoustic oscillation peak length 
scale and supernova Type la apparent magnitude data favors a spatially-flat cosmological model 
currently dominated by a time-independent cosmological constant but does not exclude time- varying 
dark energy. 
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. A number of lines of observational evidence support a "standard" model of cosmology with energy budget dominated 
by far by dark energy. Dark energy is most simply characterized as a negative-pressure substance that powers the 
observed accelerated cosmological expansion. It can evolve slowly in time and vary weakly in space, although current 
data are consistent with it being Einstein's cosmological constant. On the other hand, some argue that the observed 
\ accelerated expansion should instead be viewed as an indication that general relativity does not accurately describe 
gravitational physics on large cosmological length scales. For recent reviews see 2-8]. In what follows we assume that 
i— general relativity provides an accurate description of gravitation on cosmological length scales. 

There are many dark energy models under discussion. For recent discussions see [9l4l7|. and references therein. 
Perhaps the most economical — and the current "standard" model — is the ACDM model [18| , where the accelerated 
cosmological expansion is powered by Einstein's cosmological constant, A, a spatially homogeneous fluid with equation 
\Q • of state parameter wa = Pa/ PA = —1 (with p\ and pa being the fluid pressure and energy density). In this model the 
\Q \ cosmological energy budget is dominated by pa, with cold dark matter (CDM) being the second largest contributor. 

■ The ACDM model provides a reasonable fit to most observational constraints, although the "standard" CDM structure 
ly-^ ' formation model might be in some observational trouble (see, e.g., flfjl |20|). 

The ACDM model has a few apparent puzzles. Prominent among these is that the needed A energy density scale 
t— I ' is of order an meV, very small compared to the higher (cutoff) value suggested by a perhaps naive application of 
quantum mechanics. Another puzzle is that we happen to be observing at a time not very different from when the A 
energy density started dominating the cosmological energy budget (this is the "coincidence" puzzle) . 

If the dark energy density — that responsible for powering the accelerated cosmological expansion — slowly de- 
creased in time (rather than remaining constant like pa), the energy densities of dark energy and nonrelativistic 
matter (CDM and baryons) would remain comparable for a longer period of time, and so alleviate the coincidence 
puzzle. Also, a slowly decreasing dark energy density, that is based on more fundamental physics at a higher energy 
density scale much larger than an meV, would result in a current dark energy density scale of an meV through gradual 
decrease over the long lifetime of the Universe. Thus a slowly decreasing dark energy density could resolve some of 
the puzzles of the ACDM model [H|. 

The XCDM parametrization is often used to describe a slowly decreasing dark energy density. In this parametriza- 
tion the dark energy is modeled as a spatially homogenous (X) fluid with an equation of state parameter wx — Px/px, 
where wx < —1/3 is an arbitrary constant and px and px are the pressure and energy density of the X-fluid. When 
Wx — —1, the XCDM parametrization reduces to the ACDM model, which is a complete and consistent model. 
For any other value of wx (< — 1/3) the XCDM parametrization is incomplete as it cannot describe spatial inhomo- 
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geneities (see, e.g., |23[). For computational simplicity we study the XCDM parametrization only in the spatially-flat 
cosmological case. 

If the dark energy density evolves in time, physics demands that it also be spatially inhomogeneous. The 0CDM 
model — in which dark energy is modeled as a scalar held cf> with a gradually decreasing (in (jj) potential energy density 
V{4>) — is the simplest complete and consistent model of a slowly decreasing (in time) dark energy density. Here 
we focus on an inverse power-law potential energy density V(cf>) oc _a , where a is a nonnegative constant [2lL l24j|. 
When a — the 0CDM model reduces to the corresponding ACDM case. Here we only consider the spatially-flat 
0CDM cosmological model. 

It has been known for some time now that a spatially-flat ACDM model with current energy budget dominated 
by a constant A is largely consistent with most observational constraints (see, e.g., [25r 28]). Supernovae Type la 
(SNela) apparent magnitude measurements (e.g., (22l . [29l - [3l| ). in conjunction with cosmic microwave background 



(CMB) anisotropy data (e.g., |32j-|36|) combined with low estimates of the cosmological mass density (e.g., |37j), as 
well as baryon acoustic oscillation (BAO) peak length scale estimates (e.g., [38H4lj ) strongly suggest that we live in a 
spatially-flat ACDM model with nonrelativistic matter contributing a little less than 30 % of the current cosmological 
energy budget, with the remaining slightly more than 70 % contributed by a cosmological constant. These three sets 
of data carry by far the most weight when determining constraints on models and cosmological parameters. 

Future data from space missions will tighten the constraints (see, e.g., [42|-[44j]). However, at present, it is of great 
importance to consider independent constraints that can be derived from other presently available data sets. While 
these data do not yet carry as much statistical weight as the SNela, CMB and BAO data, they potentially can reassure 
us (if they provide constraints consistent with those from the better known data), or if the two sets of constraints 
are inconsistent this might lead to the discovery of hidden systematic errors or rule out the cosmological model under 
consideration. 

Other data that have been used to constrain cosmological parameters include galaxy cluster gas mass fraction 
(e.g., [H, |H, HH), gamma-ray burst luminosity distance (e.g. . l47H53l). larg e-scale structure (e.g., [5lT[53ln. strong 
gravitational lensing (e.g., [5J-[57j]), and lookback time (e.g., |5Ct l58l. 159. l6l|) or Hubble parameter (e.g., j62[, (63[, 
[64| . [65j]) data. While the constraints from these data are much less restrictive than those derived from the SNela, 
CMB and BAO data, both types of data result in largely compatible constraints that generally support a currently 
accelerating cosmological expansion. This gives us confidence that the broad outlines of the "standard" cosmological 
model are now in place. 

Angular size data have also been used to constrain cosmological parameters (see, e.g., [66 70]). In this paper we 
use the Bonamente et al. hereafter B06) galaxy cluster angular size versus redshift data to derive cosmological 
constraints. These measurements were determined from radio and X-ray observations. They have previously been 
used to constrain some cosmological parameters and to test the distance duality relationship of metric gravity models 
(see, e.g., 0^1)- 

In this paper we use the B06 angular size versus redshift data to constrain cosmological models not previously 
considered, and to constrain other cosmological parameters in models previously considered. We show that these 
constraints are compatible with those derived using other data. We also do a joint analysis of this angular size data 
and SNela and BAO measurements and show that including the angular size data in the mix affects the constraints, 
although not greatly so as the angular size data do not yet have sufficient weight. 

Our paper is organized as follows. In Sec. [TT] we present the basic equations of the three dark energy models we 
study. Constraints from the B06 angular diameter distances of galaxy clusters are derived in Sec. IIIII In Sec. HVI thc 
BAO data and the SNela measurements are used to constrain the dark energy models. In Sec|V]we determine joint 
constraints on the dark energy parameters from different combinations of data sets. Finally, we summarize our main 
conclusions in Sec. ED 



II. BASIC EQUATIONS OF THE DARK ENERGY MODELS 

The Friedmann equation of the ACDM model with spatial curvature can be written as 

e 2 {z- p) = n m0 (i + zf + n A + (i - n m0 - r> A )(i + zf, (i) 

where z is the redshift, E(z) = H(z)/Hq is the dimensionless Hubble parameter where Ho is the Hubble constant, and 
the model-parameter set is p = (£l m o,QA_) where f2 m o is the nonrelativistic (baryonic and cold dark) matter density 
parameter and Da that of the cosmological constant. Throughout, the subscript denotes the value of a quantity 
today. In this paper, the subscripts A, X and <f> represent the corresponding quantities of the dark energy component 
in the ACDM, XCDM and </>CDM models. 
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In this work, for computational simplicity, the spatial curvature is set to zero in the XCDM and 0CDM cases. Then 
the Friedmann equation for the XCDM parametrization is 

E 2 (z; p) = n m0 (l + z f + (1 - 0*0X1 + z f {1+Wx \ (2) 

where the model-parameter set is p = (Ono, w x)- 

In the 0CDM model, the inverse power law potential energy density of the scalar field adopted in this paper is 
V{4>) — Km.p<j)~ a , where m p is the Planck mass, and a and n are non- negative constants [13] • I n the spatially-flat 
case the Friedmann equation of the </>CDM model is 

H 2 (z;p) = (Pm + P<t>), (3) 

where H(z) = a/a is the Hubble parameter, and a(t) is the cosmological scale factor and an overdot denotes a time 
derivative. The energy densities of the matter and the scalar field are 

Prn = -^a 3 , (4) 



and 
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/>0=3^(4> 2 + «rn>- Q ), (5) 



respectively. According to the definition of the dimensionless density parameter, one has 



Q ( z \ = 8lT P m = P™ ( 6 ) 
3m%H 2 p m + p<f,' 



The scalar field <f> obeys the differential equation 

+ 3^-^m^- (a+1) =O. (7) 
a 2 * 

Using Eqs. ([3]) and (O, as well as the initial conditions described in [23 |. one can numerically compute the Hubble 
parameter H(z). In this case, the model-parameter set is p = (O m oj 

To use observational data to constrain cosmological models, we need various distance expressions. The coordinate 
distance is 



aoH y/\rik\ 



=smn 



dz' 

W) 



(8) 



where is the spatial curvature density parameter and c is the speed of light, and 



• 1 ArTl ^ f sin(V-0 x)/-y/-0 (O < 0), 
3mn \ Sit \ x) \ 

X ' ' - ] x (0 = 0), (9) 



lsinh(VOx)/^OT (O>0). 
The luminosity distance cLl and the angular diameter distance cZa are simply related to the coordinate distance as 

d L = (a r)(l+z), (10) 

and 

d A = (a r)/(l + z). (11) 



III. CONSTRAINTS FROM THE ANGULAR SIZE DATA 



X-ray observations of the intracluster medium combined with radio observations of the galaxy cluster Sunyaev- 
Zeldovich effect allow for an estimate of the angular diameter distance (ADD) d A of galaxy clusters. Here we use 
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TABLE I: Angular diameter distances of galaxy clusters from B06. 



the 38 ADDs of B06 to constrain cosmological parameters. These data can be found in Tables 1 and 2 of B06. For 
convenience, we re-collect them in Table [IJ 

There are three sources of uncertainty in the measurement of (La'- the cluster modeling error cr mo( j; the statistical 
error <7 s t a t; and, the systematic error cr sys . The modeling errors are shown in Table U and the statistical and systematic 
errors are presented in Table 3 of B06. In our analysis here we combine these errors in quadrature. Thus, the total 
uncertainty a to t satisfies of ot = cr^ od + cr s 2 tat + cr s 2 ys . 

We constrain cosmological parameters by minimizing xaddi 

2 , ^ [d%{zi;H ,p)-d>f(zi)]* 

xadd{ h o,p) = / J 2 • y 12 > 

i=X *^tot,« 

Here z% is the redshift of the observed galaxy cluster, is the predicted value of the ADD in the cosmological model 
under consideration and d° A hs is the measured value. From Xadd(^o, p) we compute the likelihood function L{Hq, p). 
We then treat Hq as a nuisance parameter and marginalize over it using a gaussian prior with Hq = 68 ± 3.5 km s _1 
Mpc -1 to get a likelihood function L(p) that is a function of only the cosmological parameters of interest. The 
best-fit parameter values p* are those that maximize the likelihood function and the 1, 2 and 3 a constraint contours 
are the set of cosmological parameters (centered on p*) that enclose 68.27, 95.45 and 99.73 %, respectively, of the 
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FIG. 1: 1, 2, and 3 a constraint contours for the ACDM model from the ADD data. The dashed diagonal line corresponds to 
spatially-flat models and the shaded area in the upper left-hand corner is the region for which there is no big bang. The star 
marks the best-fit pair (fimC^A) = (0.19, —0.62) with Xmm = 30.1. 



probability under the likelihood function. 

Figures [T|— [3] show the constraints from the ADD data on the three dark energy models we consider. Comparing 
these results to those shown in Figs. 1 — 3 of [§8[ , derived using the compact radio source angular size data of Gurvits 
et al. 66], and to Figs. 1 — 2 of (69|, derived using the FRIIb radio galaxy angular size data from Guerra et al. [67| . 
we see that the B06 galaxy cluster angular size data result in approximately comparable constraints on cosmological 
parameters as those from the two earlier angular size data sets. These ADD constraints are comparable to those from 
gamma-ray burst data ([50], Figs. 1 — 3 and 7 — 9) as well as those from Hubble parameter measurements (jfi^l, Figs. 
1-3). 

Clearly, current ADD data constraints are not very restrictive, although it is encouraging that the ADD con- 
straints on these dark energy models do not disfavor the regions of parameter space that are favored by other data. 
More importantly, we anticipate that near future ADD data will provide significantly more restrictive constraints on 
cosmological parameters. 



IV. CONSTRAINTS FROM BAO AND SNEIA DATA 



The BAO peak length scale can be used as a standard ruler to constrain cosmological parameters. Here we use 
the recent Percival et al. [78[ BAO data to constrain parameters of the ACDM and <^>CDM models and the XCDM 
parametrization. 

Percival et al. (2010) measure the position of the BAO peak from the SDSS DR7 and 2dFGRS data, determining 
r s{ z d)l Dy(z — 0.275) = 0.1390 ± 0.0037, where r s (zd) is the comoving sound horizon at the baryon drag epoch, and 
D v (z) = [(1 + z) 2 d 2 A cz/H{z)] 1 / z . By using fl m0 h 2 = 0.1326 ± 0.0063 and fl b0 h 2 = 0.02273 ± 0.00061 (here n b0 is the 
current value of the baryonic mass density parameter and h is the Hubble constant in units of 100 km s _1 Mpc ) 
from WMAP5 [35], one can further get 

/ O 7,2 \ -0134 /0 , 2 \ -0-255 

D y (0.275) = (1104 ±30)( 2f$°2_ ) Mpc, (13) 

v ' y ' V0.022737 V 0-1326 / 

as is shown in Eq. (13) of (78[. The error on Vl^h 2 is ignored in this work, as the WMAP5 data constrains Slboh 2 to 
0.5 %. 

Our results for the ACDM model and the XCDM parametrization agree very well with the Ref. [78[ results shown in 
their Fig. 5. Our results for the 0CDM model are shown in Fig.UJ BAO data primarily constrains fl m o significantly, 
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FIG. 2: 1, 2, and 3 <j constraint contours for the XCDM parametrization from the ADD data. The dashed horizontal line 
at ujx = — f corresponds to spatially-flat ACDM models. The star marks the best-fit pair (QmOjWx) = (O.Of , — 0.12) with 
xLn = 30.2. 




FIG. 3: 1, 2, and 3 a constraint contours for the 0CDM model from the ADD data. The horizontal axis at a — corresponds 
to spatially-flat ACDM models. The star marks the best-fit pair (O m o, a) = (0.54, 5) with Xmln = 37.3. 

leaving 17a, wx and a almost unconstrained (see, e.g., [79j]). The results from the BAO data [78| are most directly 
comparable to those derived from the earlier BAO data of Ref. jsjj. Comparing to the constraints shown in Figs. 
2 — 4 of Ref. [7^], one sees that the Percival et al. (2010) data results in slightly more restrictive constraints than the 
Eisenstein et al. (2005) data. 

Type la supernovae are standardizable candles that can be used to constrain cosmological parameters. Here we use 
the recent Union2 compilation of 557 SNela (covering a redshift range 0.015 < z < 1.4) of Amanullah et al. [81J to 
constrain parameters of the ACDM and 0CDM models and the XCDM parametrization. 
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mO 

FIG. 4: 1, 2, and 3 cr constraint contours for the </>CDM model from the BAO data. The horizontal axis at a = corresponds 
to spatially-flat ACDM models. The star marks the best-fit pair (f2 m o, a) = (0.32, 2.01) with xLn = 0.169. 

Cosmological constraints from SNela data are obtained by using the distance modulus fJ,(z). The theoretical 
(predicted) distance modulus is 

/Zth(z;Pi/io) = 51og 10 [D L (z;p)] +/x 0) (14) 
where (1q — 42.38 — 5 log 10 h and the Hubble-free luminosity distance is given by 

D L (z;p) = ^d L = (l + z) f^T-T- (15) 
c J E(z';p) 

The best-fit values of cosmological model parameters can be determined by minimizing the % 2 function 

557 

2 / \ Mth.iO^P'Mo) - Mobs,i( z «) 

Xsw(P> Mo) = >. 5 , (16) 

— err, 

where Mobs iC^i) is the distance modulus obtained from observations and er^ is the total uncertainty of the SNela data. 
The zero-point /io is treated as a nuisance parameter and marginalized over analytically [82T - l84l | . The covariance matrix 
with or without systematic errors can be found on the web (http:/ /supernova. lbl.gov/Union). 

Our results for the ACDM model and the XCDM parametrization agree very well with the Ref. (si[ results shown in 
their Figs. 10 and 11. The constraints on <^CDM model parameters from these data are shown in Fig. [5] Comparing 
to Fig. 4 of Ref. [4(|, we see that the constraints from the Ref. [8l[ data with systematics errors are approximately 
comparable to those from the earlier Ref. [85[ data for 307 SNela without systematic errors. Like the case for the 
BAO data, the SNela data constraints are also fairly one dimensional, tightly constraining one combination of the 
cosmological parameters while only weakly constraining the "orthogonal" combination. 



V. JOINT CONSTRAINTS 



Figures [B]— [5] show the constraints on the cosmological parameters for the ACDM and </>CDM models and the XCDM 
parametrization from a joint analysis of the BAO and SNela data, as well as from a joint analysis of the BAO, SNela 
and ADD data. With the inclusion of systematic errors in the analysis of the SNela data of Ref. 
BAO and SNela constraints (thin solid contours in Figs. [B] 

Ref. 50] that made use of the smaller SNela data set of Ref. 85] that did not include systematic errors. 



the new joint 

are similar to the earlier ones shown in Figs. 4 — 6 of 
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FIG. 5: 1, 2, and 3 cr constraint contours for the </>CDM model from the SNela data. The horizontal axis at a = corresponds 
to spatially-flat ACDM models. Thin solid lines (best fit at Q. m o = 0.27 and a — 0.0 with Xmm = 543, marked by "x" ) exclude 
systematic errors, while thick solid lines (best fit at Q m o = 0.27 and a — 0.0 with Xmin — 531, marked by "0") account for 
systematics. 



Model 


BAO + SNela 


ADD + BAO + SNela 


ACDM 


0.24 < fi m0 < 0.33 
0.5 < n A < 0.97 


0.24 < fi m0 < 0.33 
0.46 <fl A < 0.93 


XCDM 


0.24 < n m o < 0.33 
-1.30 < u) X < -0.80 


0.24 < n m0 < 0.33 
-1.25 < u)x < -0.77 


<^>CDM 


0.24 < Q m o < 0.33 
< a < 0.73 


0.24 < Q m0 < 0.33 
0.01 < a < 0.89 



TABLE II: Two standard deviation bounds on cosmological parameters. 

Figures 191— 1111 display the one-dimensional marginalized distribution probabilities of the cosmological parameters 
for the three cosmological models considered in this work, derived from a joint analysis of the BAO and SNela data, 
as well as from a joint analysis of the BAO, SNela and ADD data. The marginalized 2 a intervals of the cosmological 
parameters are presented in Table [XXJ 

The combination of BAO and SNela data gives tight constraints on the cosmological parameters. Adding the 
currently-available galaxy cluster ADD data to the mix does shift the constraint contours, however the effect is not 
large. Current ADD data do not have enough weight to significantly affect the combined BAO and SNela results. 
The ADD data have approximately the same weight as currently-available gamma-ray burst luminosity measurements 
(HH, Figs. 4— 6 and 10— 12). 

Clearly, the observational data considered here are very consistent with the predictions of a spatially-flat cosmo- 
logical model with energy budget dominated by a time-independent cosmological constant. However, the data do not 
rule out time-evolving dark energy, although they do require that it not vary rapidly. 

VI. CONCLUSION 

We have shown that the galaxy cluster angular size versus redshift data of Ref. |lj can also be used to constrain 
dark energy model cosmological parameters. The resulting constraints are compatible with those derived from other 
current data, thus strengthening support for the current "standard" cosmological model. The ADD constraints are 
approximately as restrictive as those that follow from currently available gamma-ray burst luminosity data, strong 
gravitational lensing measurements, or lookback time (or Hubble parameter) observations. They are, however, much 
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FIG. 6: Thick (thin) solid lines are 1, 2, and 3 a constraint contours for the ACDM model from a joint analysis of the BAO and 
SNela (with systematic errors) data, with (and without) the ADD data. The cross ("+") marks the best-fit point determined 
from the joint sample without the ADD data at f2 m o = 0.28 and Qa = 0.76 with Xmin = 531. The star ("*") marks the best-fit 
point determined from the joint sample with the ADD data at Q m o = 0.28 and Qa = 0.72 with Xmin = 565. The dashed sloping 
line corresponds to spatially-flat models. 

less restrictive than those that follow from a combined analysis of BAO peak length scale and SNela apparent 
magnitude data. 

The spatially- flat ACDM model, currently dominated by a constant cosmological constant, provides a good fit to 
the data we have studied here. However, these data do not rule out a time-evolving dark energy. 

We anticipate that near-future angular size data will provide significantly more restrictive constraints on cosmolog- 
ical parameters. In conjunction with other observations, this angular size data will prove very useful in pinning down 
parameter values of the "standard" cosmological model. 
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FIG. 10: One-dimensional marginalized distribution probabilities of the cosmological parameters for the XCDM parametriza- 
tion. Thick (thin) lines are results from a joint analysis of the BAO and SNela (with systematic errors) data, with (and without) 
the ADD data. 
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FIG. 11: One-dimensional marginalized distribution probabilities of the cosmological parameters for the (j>CDM model. Thick 
(thin) lines are results from a joint analysis of the BAO and SNela (with systematic errors) data, with (and without) the ADD 
data. 



